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Abstract 



In a general two-Higgs-doublet model (2HDM), without the ad hoc discrete symme- 
l/"") . tries to prevent tree-level flavor-changing-neutral currents, an extra phase angle in the 

charged-Higgs-fermion coupling is allowed. We show that the charged-Higgs amplitude 
interferes destructively or constructively with the standard model amplitude depend- 
ing crucially on this phase angle. The popular model I and II are special cases of our 
analysis. As a result of this phase angle the severe constraint on the charged-Higgs 
boson mass imposed by the inclusive rate of b — > 57 from CLEO can be relaxed. We 
also examine the effects of this phase angle on the neutron electric dipole moment. 
Furthermore, we also discuss other constraints on the charged-Higgs-fermion couplings 
qj' coming from measurements of B° — B° mixing, po, and Rf,. 

^" 

I. Introduction 

X 

One of the most popular extensions of the standard model (SM) is the two-Higgs-doublet 
model (2HDM) [Jj , which has two complex Higgs doublets instead of only one in the SM. The 
2HDM allows flavor- changing neutral currents (FCNC), which can be avoided by imposing 
an ad hoc discrete symmetry M . One possibility to avoid the FCNC is to couple the fermions 
only to one of the two Higgs doublet, which is often known as model I. Another possibility 
is to couple the first Higgs doublet to the down-type quarks while the second Higgs doublet 
to the up-type quarks, which is known as model II. This model II has been very popular 
because it is the building block of the minimal supersymmetric standard model. The physical 
content of the Higgs sector includes a pair of CP-even neutral Higgs bosons H° and h , a 
CP-odd neutral boson A, and a pair of charged-Higgs bosons H . 

Models I and II have been extensively studied in literature and tested experimentally. 
One of the most stringent tests is the radiative decay of B mesons, specifically, the inclusive 
decay rate of b — > 57, which has the least hadronic uncertainties. The SM rate of b — > 57 
including the improved leading-order logarithmic QCD corrections is predicted H to be 
(2.8±0.8) x 10~ 4 , of which the uncertainty mainly comes from the factorization scale and from 



the next-to-leading order corrections. [] In the 2HDM, the rate of b — > s 7 can be enhanced 
substantially for large regions in the parameter space of the mass M^± of the charged-Higgs 
boson and tan (3 = v 2 /vi, where V\ and v 2 are the vacuum expectation values of the two Higgs 
doublets. CLEO published a result of b — > s 7 inclusive rate of (2.32 ± 0.57 ± 0.35) x 10~ 4 in 
1995 §, which is recently updated to (3.15 ±0.35 (stat) ±0.32 (sys) ±0.26 (mod)) x 10" 4 in 
1998 @. ALEPH also published a result of (3.11 ± 0.80 (stat) ± 0.72 (sys) ) x 10" 4 §. The 
95%CL limit published by CLEO is also updated to 2 x 10~ 4 < B(b -> s 7 ) < 4.5 x 10" 4 [§. 
The data is now more consistent with the SM prediction than before. Hence, the experimental 
result puts a rather stringent constraint on the charged-Higgs boson mass M H ± and tan/3. 
In model II, the constraint is M#± > 350 GeV for tan/3 larger than 1, and even stronger for 
smaller tan/? ||. 

Recently, there have been some studies |9|, |T(J on a more general 2HDM without the 
discrete symmetries as in models I and II. It is often referred as model III. FCNC's in 
general exist in model III. However, the FCNC's involving the first two generations are 
highly suppressed from low-energy experiments, and those involving the third generation is 
not as severely suppressed as the first two generations. It implies that model III should be 
parameterized in a way to suppress the tree-level FCNC couplings of the first two generations 
while the tree-level FCNC couplings involving the third generation can be made nonzero as 
long as they do not violate any existing experimental data, e.g., B° — B° mixing. 

In this work, we simply assume all tree-level FCNC couplings to be negligible. Even 
though in such a simple model the couplings involving Higgs bosons and fermions can have 
complex phases e %e . The effects of such extra phases in b — ► s 7 have been noticed in Ref . [[iT| . 
In this paper, we shall study carefully the constraint on the phase angle in the product, 
AttAftfe, of Higgs-fermion couplings (see below) versus the mass of the charged-Higgs boson 
from the CLEO data of b — > 57. We shall show that in the calculation of b — > s 7 the 
charged-Higgs amplitude interferes destructively or constructively with the SM amplitude 
depending crucially on this phase angle and less on the charged-Higgs mass. The usual model 
I and II are special cases in our study. We shall also show that the previous constraints on 
the charged-Higgs mass and tan (3 imposed by the CLEO data can be relaxed because of the 
presence of this extra phase angle. There are other processes in which the effects of the phase 
angle can be seen. One of these that we study in this report is the neutron electric dipole 
moment. In addition, we also discuss the constraints from experimental measurements of 
B° — B° mixing, p , and Rb. 

The organization is as follows. In the next section we describe the content of the general 
2HDM and write down the Feynman rules for model III. In Sec. Ill, we describe briefly the 
effective hamiltonian formulation for the decay of b — > 37 and derive the Wilson coefficients 
in model III. We present our numerical results for b — > s 7 and study the case of neutron 
electric dipole moment in Sec. IV. In Sec V, we discuss other experimental constraints from 
measurements of B° — B° mixing, po, and Rb- Finally, we conclude in Sec. VI. 



*Thc NLO order calculations for the SM and 2HDM I and II are available very recently (§. The SM 
result is (3.29 ± 0.33) x 10 -4 , which is consistent with the LO calculation. However, the NLO calculation is 
not available for 2HDM III and, therefore, we will use the LO result consistently throughout the paper. 
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II. The General Two-Higgs-Doublet model 



In a general two-Higgs-doublet model, both the doublets can couple to the up-type and 
down-type quarks. Without loss of generosity, we work in a basis such that the first doublet 
generates all the gauge-boson and fermion masses: 




i> = 



X 



where v is related to the W mass by Mw = \v. In this basis, the first doublet 4>i is the same 
as the SM doublet, while all the new Higgs fields come from the second doublet 02- They 
are written as 



, J_ / y/2G+ 




V2H+ 
X°2 + iA° 



(2) 



where G° and G are the Goldstone bosons that would be eaten away in the Higgs mechanism 
to become the longitudinal components of the weak gauge bosons. The H ± are the physical 
charged-Higgs bosons and A is the physical CP-odd neutral Higgs boson. The Xi an d X2 are 
not physical mass eigenstates but linear combinations of the CP-even neutral Higgs bosons: 



Xi 

xl 



H° cos a — h° sin a 
H° sin a + h° cos a 



(3) 
(4) 



where a is the mixing angle. In this basis, there is no couplings of x\Z% an d X2W + W . 
We can write down|TD] the Yukawa Lagrangian for model III as 



£y = VijQiL(j>iU jR + riijQiLfaD jR + &Q iL faU jR + C, O.io^n,!! + h.c. 



(5) 



icT20i,2 , Vij D an d £ij D are > m general, nondiagonal 



where i,j are generation indices, fag 
coupling matrices, and Qil is the left-handed fermion doublet and UjR and Dj R are the 
right-handed singlets. Note that these Qn, Ujr, and Dj R are weak eigenstates, which can 
be rotated into mass eigenstates. As we have mentioned above, 0i generates all the fermion 
masses and, therefore, ^V U,D wm become the up- and down-type quark-mass matrices after 
a bi-unitary transformation. After the transformation the Yukawa Lagrangian becomes 



C Y = -UM V U - DM D D - 



2M W 



(# u cos a - h" sin a) UM V U + DM D D 



+ 
+ 



2M W 
9 



G° (UMuj 5 U - DMr>j 5 D 
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V2M W G+UVCKM 
H° sin a + h° cos a 

71 



Mu\{\ 



7 



u 

D 



u 



3 



+ 



+ J D^ D |(l + 7 5 )+^ t |(l-7 5 



D 



V2 
H + 'U 



7 



U-D + 7 B)_ pKl-^) d 



*W D i(l + 7 5 ) 



S l/ CKM2V i ' 



^VckmKI-7 5 



U 



(6) 



where [/ represents the mass eigenstates of u, c, t quarks and D represents the mass eigen- 
states of d, s, b quarks. The transformations are defined by M^d = diag(m Uj d, m CjS , m tt b) = 
-^(£u,Dyv U ' D (Ku,D), i U ' D = (C UiD )^ u ' D (TZ UiD ). The Cabibbo-Kobayashi-Maskawa matrix 
is Vckm = (CuY(C D ). 

The FCNC couplings are contained in the matrices £ U ' D '. A simple ansatz for £ ,D would 
bei 



V2M W ^ 

by which the quark-mass hierarchy ensures that the FCNC within the first two generations 
are naturally suppressed by the small quark masses, while a larger freedom is allowed for 
the FCNC involving the third generations. Here Ay's are of order 0(1) and unlike previous 



studies H [UJ they can be complex, which give nontrivial consequences different from pre- 
vious analyses based on model I and II. An interesting example would be the inclusive rate 
of b — ► S7 that we shall study next. Such complex Ay's allow the charged-Higgs amplitude 
to interfere destructively or constructively with the SM amplitude. As we have mentioned, 
models I and II are special cases in our study and so the previous constraints |§ imposed 
on the charged-Higgs mass and tan/5 by the CLEO data can be relaxed by the presence 
of the extra phase angle. Other interesting phenomenology of the complex A^-'s includes 
the electric dipole moments of electrons and quarks |Tj| as a consequence of the explicit CP 
violation due to the complex phase in the charged-Higgs sector. For simplicity we choose 
£ U ' D to be diagonal to suppress all tree-level FCNC couplings and, consequently, the A^-'s 
are also diagonal but remain complex. Such a simple scenario is sufficient to demonstrate 
our claims. 



III. Inclusive B —> X s j 



The detail description of the effective hamiltonian approach can be found in Refs. M, [14 . 



Here we present the highlights that are relevant to our discussions. The effective hamiltonian 
for B — > A s 7 at a factorization scale of order 0(m&) is given by 



n 



eff 



(8) 



The operators Qi can be found in Ref.|§, of which the Qi and Q2 are the current-current 
operators and Q3 — Qq are QCD penguin operators. Q 7l and Qsg are, respectively, the 
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magnetic penguin operators specific for b — > 57 and b — > sg. Here we also neglect the mass 
of the external strange quark compared to the external bottom-quark mass. 

The factorization in Eq.(|8|) facilitates the separation of the short-distance and long- 
distance parts, of which the short-distance parts correspond to the Wilson coefficients 
and are calculable by perturbation while the long-distance parts correspond to the operator 
matrix elements. The physical quantities based on Eq. (H) should be independent of the 
factorization scale \l. The natural scale for factorization is of order m b for the decay B — > 
X s 7. The calculation of the Cj(//)'s divides into two separate steps. First, at the electroweak 
scale, say M w , the full theory is matched onto the effective theory and the coefficients 
Ci(M\y) at the H^-mass scale are extracted in the matching process. In a while, we shall 
present these coefficients Ci(M w ) in our model III. Second, the coefficients Ci(M w ) at the W- 
mass scale are evolved down to the bottom-mass scale using renormalization group equations. 
Since the operators Qj's are all mixed under renormalization, the renormalization group 
equations for Cj's are a set of coupled equations: 

C(fi) = U(ti,M w )C(M w ) , (9) 

where U(ji,Mw) is the evolution matrix and C(/jl) is the vector consisting of Cj(ju)'s. The 
calculation of the entries of the evolution matrix U is nontrivial but it has been written down 
completely in the leading order B. The coefficients Ci(fi) at the scale 0{m b ) are given by 



Cj(jj) = E fc /' (j = l,-,6) , (10) 
1=1 

CV 7 (/i) = V %C 7 ,(M W ) + I( V % - V %)C 8 G(M w )+C 2 (M w )j2hiV ai , (11) 



i=i 

CM = v^C SG (M w ) + C 2 (M w )j2h i V a ' , (12) 

i=i 

with r] = a s (M w )/a s (fi). The aj's, fc^'s, h^s, and hi's can be found in Ref. ||. 

Once we have all the Wilson coefficients at the scale 0(m b ) we can then compute the 
decay rate of B — > X s / y. The decay amplitude for B — ► X s ^ is given by 

A(B -> X sl ) = -9Zv; s V tb C 7 M(Q 7 ,) , (13) 

in which we use the spectator approximation to evaluate the matrix element (Q 7l ) and 
rriB — m b . The decay rate of B — > X s 7 is given by 

Gl\V * s V tb \ 2 a Prn ml 
32vr 4 

where C 7j (m b ) is given in Eq. (TP]). Since this decay rate depends on the fifth power of m&, 
a small uncertainty in the choice of mj will create a large uncertainty in the decay rate, 
therefore, the decay rate of B — > X s 7 is often normalized to the experimental semileptonic 
decay rate as 



T(B - X s7 ) = " F|Vf C" emmfe |C 77 (m b )| 2 , (14) 



T(S -> X c ez/ e ) (V^bl 2 7rf(m c /m b ) 



C 7j (m b )\\ (15) 
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where f(z) = 1 - 8z 2 + 8z 6 - z 8 - 2Az 4 In z. 

The remaining task is the calculation of the Wilson coefficients Cj(M^) at the 1^-mass 
scale. The necessary Feynman rules can be obtained from the Lagrangian in Eq. (|6|). As 
we have mentioned, we assume all tree-level FCNC couplings negligible and, therefore, the 
neutral-Higgs bosons do not contribute at tree level or at one-loop level. The only contri- 
butions at one-loop level come from the charged-Higgs bosons H , the charged Goldstone 
bosons G ± , and the SM W ± bosons. 

The coefficients Cj(M^) at the leading order in model III are given by 



Cj{M w ) 
C 2 (M W ) 

C 7j (M w ) 
C 8G (M W ) 




1 



(j 



Ajxt) 
2 

D(xt) 



1,3,4,5,6) 

My) 



6 

D(y) 



X tt \ 2 + B(y)X tt X 



66 



X tt \ 2 + E(y)X tt X 



2 6 

where x t = m 2 /M^, and y = m 2 /Mf I ±. The Inami-Lim functions [jl 

"8x 2 + 5x-7 (3s 2 - 2x) lux 



bb 



(16) 
(17) 



(19) 



are given by 



A(x) 
B(y) 

D(x) 

E(y) 



X 



y 



X 



y 



12(x - l) 3 
5y-3 



(3y 



2(x - l) 4 
-2) In 



x 



5x-2 



4(x - l) 3 

y-3 



+ 



6(y-l) 3 
3x lnx 



+ 



2(x- l) 4 
In y 



4(y-l) 2 2(3/ -1 



(20) 
(21) 
(22) 
(23) 



The SM results for the Wilson coefficients Ci(M w ) for i = 1, 6 are the same as in Eqs. ( jig ) 
and (|ITD, while C 7l (M w ) and C 8 g(M w ) only have the first term as in Eqs. ( p~8|) and ([T9|), 
respectively. Thus, we already have all the necessary pieces to compute the decay rate of 
B -> X s7 . 

Before we leave this section we would like to emphasize that the expressions for Ci(Mw) 
in Eqs. ( |T6"D - ( |I~9"D obtained for model III can be reduced to the results of models I and II 
by the following substitutions: 



and 



X u — > cot (3 and 



X u — > cot (3 and X bb 



-> cot (3 (for model I), 
— taxi (3 (for model II). 



(24) 
(25) 



IV. Results 



We use the following inputs 0, [16], [T?| for our calculation: 



GcV. 



KM 



v, 



cb\ 



1/133, and a s (M z ) 



m t = 173.8 GeV, M w = 80.388 
= 0.95, m c /m b = 0.3, and B(b -> ce^z/) = 10.45 ± 0.21%, a om (m b ) ~ 
0.119 and a 1-loop a s is employed. The branching ratio B(B — > X s 7) 
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is calculated using Eq. (p~5|) . The free parameters are then M H ±, X u , and A^, as in Eqs. ( |l"8D 
and (fT - 



Since the term proportional to XtAbb is, in general, complex we let Xu^bb = \KtXbb\e ■ 
We show the contours of the branching ratio in the plane of 9 and M#± for | At* A&b | = 3, 1, 0.5 
in Fig. [1] (a), (b), and (c), respectively. The contours are symmetric about 9 = 180°. The 
contours are B = (2, 2.8, 4.5) x 10~ 4 , which correspond to 95%CL lower limit, the SM value, 
and the 95%CL upper limit. The value of |A&&| is set at 50 as preferred in the Rb constraint 
that will be shown in the next section. The corresponding values of \X U \ are 0.06,0.02, and 
0.01, which satisfy the constraint from the B° — B° mixing, as will also be discussed in the 
next section. Here the term proportional to |A^| 2 is not crucial because the coefficient of 



\X U \ 2 is small compared with other two terms in Eqs. (|18|) and (19). 

The results of the conventional model II (which can be obtained from our general results 
by the substitution: X tt — > cot/?, A^ — > — tan/?) can be read off from Fig. |l|(b) at 9 = 180°. 
The b — > S7 data severely constrains M H ± > 350 GeV at 95%CL level, because at 9 = 180° 
the SM amplitude interferes entirely constructively with the charged Higgs-boson amplitude. 
It is obvious that at other angles the mass of the charged Higgs-boson mass is less constrained, 
especially, in the range 9 = 50° — 90° the entire range of charged Higgs-boson mass is allowed 
by the b — ► 57 constraint as long as |A tt A^| < 1. However, when | A it A^t, | is getting larger, 
say 3, (see Fig. |l](a)) the allowed range of charged Higgs-boson mass becomes narrow. 
This is because the charged Higgs-boson amplitude becomes too large compared with the 
SM amplitude. On the other hand, when |AtiA&&| becomes small the allowed range charged 
Higgs-boson mass is enlarged, as shown in Fig. |l|(c). The significance of the phase angle 9 
is that the constraints previously on M H ± and tan/3 are evolved into 9, M H ±, X u , and Xbb, 
where we do not need to impose \X tt \ = l/\\bb\, as in model II. The previous tight constraint 
on mn± is now relaxed down to virtually the direct search limit of almost 60 GeV at LEPII 
0- 

The phase 9 of X u Xbb can give rise to the neutron electric dipole moment (NEDM). The 
physics involved can be understood as follows. First, at the electroweak scale the phase 9 
induces the CP violating color dipole moment (CDM) of the b quark. Second, the CDM of b 



quark evolves by renormalization to the scale at m& and turns into the Weinberg operator 1 19 



i.e. the gluonic CDMj2(|) when the 6-quark field is integrated away. Finally, this gives 



NEDM at the nucleon mass scale: 

d g n = g 3 s (n)C g (fi)(O g (n)), where O g = \f abc e & ^G^G\ 5 Gl u . (26) 

Weinberg suggested the hadronic scale fi to be set at the value such that g s (^) = 4:Tt/\^6. 
Instead we choose fi at the nucleon mass. The hadronic matrix element (O g (fi)) is very 
uncertain. A typical estimate from the naive dimension analysis (NDA)[|2TJ relates the 
matrix element to the chiral symmetry breaking scale M x = 2ttF 7T = 1.19 GeV, 

O g = ■ (27) 

The parameter £ g is set to be 1 in NDA, but other calculations result in different £ g . QCD sum 
rule performed by Chemtob[^] gives £ g = 0.07. Scaling argument by Bigi and UraltsevpS 
yields a value £ 5 = 0.03. We choose £ 5 = 0.1 for our analysis. The Wilson coefficient of the 
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Weinberg operator C„ evolves according to the RG equation p4] and matches [53, E(J that 



induced by the CDM Ct, of the b quark at the scale m^. Our definitions of Wilson coefficients 



follow the notation in Ref. 126 



54 . .54 



CM = 4- 2 C b ( mi ) { ^) " f ^fi " . (28) 
32tH \oc s {m c ) J \ a s (p) J 

The CDM of the b quark comes from the CP violation of the charged Higgs coupling at the 
electroweak scale and at the scale m;, it is given by 

V2G Ft ^ x ,2 
16vr 2 




C b(m b ) = Im(A tt A b6 )~# [ ] { ) . (29) 



where the function H is 



2(1-2/) 

Note that H(l) = 1 when M#± = m t . Numerically. 



„ M = ly( y - 3 - 2 ^y] . ,30) 




< = 10- 25 e.cm Im(A B A M ) ( ] I 7?T I » ( TtM ■ « 1 1 

The experimental limit, 

d n < l(r 25 e-cm , (32) 

places an upper bound |Im(A it Afefe)| ^ 1 on the coupling product for our choice of parameters, 
£ ff = 0.1, fi = m n when Mh± ~ m t . The bound is sensitive to uncertainties in [i and £ s , but 
not much in M H ±. The function value H decreases only by a factor of 1.6 as the charged 
Higgs mass varies from 50 GeV to 200 GeV. 

In Fig. 1, the constraint on the M#± versus 9 is given by the shaded areas which are 
excluded by the NEDM measurement. 

For the case of rather large lA^A^I ^> 1, the phase becomes restricted to the forward 
region 9 ~ or the backward region ir. However, the backward region {9 ~ tt) is not 
preferable for Mjj± < 500 GeV due to the constraint from b — > sj. If the charged Higgs 
boson is this light with large couplings to the b and t quarks, the NEDM analysis requires 
a small phase in the forward region. On the other hand, when jA^A^I < 0.7, the NEDM 
constraint becomes ineffective and the constraint from b — > s'j remains useful. 

Other places to look for the effects of this angle 9 include other b — > s, d decays, CP 



violation effects in b — ► 57 [O], b — > s££, and the electric dipole moments of fermions via a 



2- loop mechanism [13 



On the other hand, this phase angle 9 will not show up in other existing constraints like 
Po, Rf,, and flavor-mixing. The previous argument that the 2HDM only has a very narrow 
window left to accommodate all the constraints from B(b — > 57), po, Rf,, and flavor- mixing 
is now not true because of the possible phase angle in model III that we are considering. 
The narrow window on M#± opens up. We shall summarize the other constraints on X tt , 
Xbb, and Higgs masses in the next section. 



8 



V. Other Constraints 



Direct searches for Higgs bosons in 2HDM at LEPII |1| place the following limits on Higgs 
boson masses: 

M h o > 77 GeV , M A > 78 GeV , M H ± > 56 - 59 GeV , (33) 

where the M^o and Ma mass limits are obtained by combining the four LEP experiments 
but no combined limit on M B ± is available [|18|]. We shall then discuss other constraints from 
precision measurements. 



V.A K° - K°, D° - D°, and B° - B° 



flavor-mixing processes can occur via tree-level, penguin, 
TOfl . One particular argument against the model III is that 



These F° — F° (F = K, D, B 
and box diagrams in model III 
it allows FCNC at the tree-level, but with a lot of freedom in picking the parameters Ay 
it certainly survives all the present FCNC constraints. The tree-level diagrams for these 
AF = 2 processes can be eliminated by choosing X U i, Xdj very small. Actually, in our study 
we have set \j = (i ^ j), therefore, all tree-level FCNC diagrams are eliminated and so 
do the penguin diagrams. However, there are important contributions coming from the box 
diagrams with the charged Higgs boson. Naively, to suppress the charged Higgs contribution 
we need to increase the charged Higgs mass or decrease \ u ■ We shall obtain a set of bounds 
using the experimental measurement Xd of B° — B° in the following ( K° — K° and D° — D° 
mixings are small in our model because of the mass hierarchy choice of in Eq. (^)). 



The quantity that parameterizes the B° — B° mixing is 



AniB 



B 



_£.IT/*I 2 

67r2 l^l 



Vtb\ f B B B m B riBTBM w (I ww + I WH + I HH ) 



(34) 



where 27 



Iww 
Iwh 
Ihh 



where x 



3 — 9x 6x 2 log x 

1 + 7 7TTT + 



xy\\u\ 2 
xy\\ tt \ 4 



x-iy (x-iy 

(Az - 1) fogy 



3 logx 



+ 



x - 4 



4 



2(l-y) 2 (l-z) 2(1 -x) 2 (l-z) 2(l-x)(l-y) 
l + y 2y\ogy 



(1 //)-' (1 //) : 



y 



m 2 /M w , 
T e use th 
0.55, and Xd 



m 2 /M 2 



z 



M w /M n:± , and the running top mass mt = m t (mt) 



166 GeV. We use these inputs P,|T|,|]: \V tb \ = 1, f 2 B B B = (0.175 GeV) 2 (1.4), m B = 5.2798 



GeV, t]b - 
from 0.004 to 0.013 
x d and it gives \y ti \ 
IT 



17 



tb\ — 

= 0.734 ± 0.035, t b = 1.56 ps. Since the allowable range of \V t d\ is 
we use a central value for \V t d\ obtained using the central value of 



0.0084 (which is the central value given in the Particle Data Book 98 
We then obtain bounds on X tt and M B ± by the 2<j limit of Xd assuming the only error 
comes from Xd measurement (see Fig. 0): 



M H ± > 77 (60) GeV for |A«| ^ 0.3 (0.28) 



(35) 
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which means virtually no limit on the charged Higgs mass if \X a \ < 0.28, because the present 
direct search limit on charged Higgs boson is about 56-59 GeV (Eq.(|33"D). We have improved 
the results in Ref. |28], [UJ because we are using an updated value of x^. In the context of 
model I and II the bound is M H ± > 77(60) GeV for tan/3 > 3.3(3.6). For tan/3 gets close to 
1, M H ± > 1 TeV. 



V.B 



Po 



p was introduced to measure the relation between the masses of W and Z bosons. In the SM 
p = Myy/Mz cos 9 W = 1 at the tree- level. However, the p parameter receives contributions 
from the SM corrections and from new physics. The deviation from the SM predictions is 
usually described by the parameter p defined by 



Po 



M 2 W 



pM% cos 2 e v 



(36) 



where the p in the denominator absorbs all the SM corrections, among which the most 
important SM correction at 1-loop level comes from the heavy top-quark: 



1 + 



3G, 



m 



t i 



(37) 



p ~ 1 + Ap t0 p 

in which Ap top is about 0.0095 for rrit = 173.8 GeV. By definition po = 1 in the SM. The 

i n nm v 

(38) 



reported value of po is |Tj| 



0.9996 ^0013 (2a) 



In terms of new physics (2HDM here) the constraint becomes: 

- 0.0017 < Ap 2H DM < 0.0013 . 



(39) 



In 2HDM po receives contribution from the Higgs bosons given by, in the context of model 

Gp 



in, m n, [kj 

Ap2HDM 

where 



8^ 



7T Z 



sin 2 aF(M H ± } M A} M H o) + cos 2 aF(M H ±,M A , M h o) 



(40) 



.F(mi, m2, 7713) = m 1 — 



' log -4 



77T.J — 77T.2 



2 2 lQ g ~\ + 2 2 lQ g "I 

mf — ITI3 \ m 3 J m 2 ~ m 3 \ m 3 



Since po is constrained to be around 1 we have to minimize the contributions of Ap 2 HDM- 
Without loss of generosity we set a = 0, which means that the heavier neutral Higgs H° 
decouples and the first Higgs doublet can be identified as the SM Higgs doublet, while the 
second Higgs doublet is the source of new physics. The leading behavior of Ap 2 HDM scale 
as Mfj ± and, therefore, the constraint of p in Eq. ( p8|) puts an upper bound on M H ±. 
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Actually, if the charged Higgs mass M H ± is between Ma and M h o the Ap 2 HDM is negative. 
However, this is not the favorite scenario because in the case of Rb the experimental result 
prefers Ma — M h o ~ 80 — 120 GeV, that will be discussed in the next subsection. In this 
case Ma — M h o, Ap 2 HDM is positive and, therefore, we want to keep it small. Using Eq. ( f40"D 
for Ma — M h o = 80 — 120 GeV, the charged Higgs mass is constrained to be 

M H ± < 180 - 220 GeV . (41) 



V.C R b 

Rb was about +3.7a above the SM value a few years ago, but now the deviation is reduced to 
+ lcr after almost all LEP data have been analyzed [16|]. i?^ xp still places a constraint on the 



2HDM, though it is much less severe than before. This is because only a narrow window exists 
in the neutral Higgs bosons that does not decrease Rb while the charged-Higgs boson always 
decreases Rb- We shall divide the discussion into two parts: neutral-Higgs contribution and 
charged-Higgs contribution. 

According to Ref. [EDI the contribution from the neutral Higgs boson is positive in a 
narrow window of 20 GeV < Ma — M^o < 120 GeV and is negative otherwise. Since 
the charged Higgs boson contribution always decreases Rb, it makes more sense to require 



the neutral Higgs contribution to be positive. Here we adapt the formulas in Ref. p9| to 
model III. First, the contribution from the neutral Higgs bosons only depends on |A&{,| and 
the masses of the neutral Higgs bosons. Again without loss of generosity, we set the scalar 
Higgs-boson mixing angle a = in order to decouple the heavier H°. We show the resultant 
Rb due to the presence of the neutral Higgs bosons in Fig. |3](a) for \\ b b\ = 30, 50, 70, where 
Rf M = 0.2158, i^ xp = 0.21656 ± 0.00074 and the la is taken to be the standard 
deviation of the experimental result. In Fig. ^(a) the horizontal lines represent the Rf M , 
+ la, and +2cr values. The R e b xp is almost at the +la line. If we allow only la value below 
i?6 Xp , we need M h o « Ma ~ 80 — 120 GeV with a fairly large |Aw,|. For |Aw,| as large as 70 
the enhancement can be as large as +la at M#± = 80 GeV. On the other hand, if we allow 
2a below R e b * v , then we can have all the range of M h o « Ma > 80 GeV, as can be seen in 
Fig. 0(a). At any rate, the preferred scenario is M h o w Ma = 80 — 120 GeV with a fairly 
large | A&& | . How large |A&&| should be? It depends on the charged Higgs contribution as well. 

Since the charged Higgs-boson contribution is always negative, we want to make it as 
small as possible. This contribution depends on |A tt |, |A^|, and M H ±. The effect on R b due 
to the presence of the charged Higgs boson is shown in Fig. |3|(b) for |A&&| = 30,50,70 and 
\X U \ = 0.05. In Fig. 0(b) the horizontal lines represent the Rf M and ±la. The R^ p is 
very close to +la line. It is clear from the graph that because we do not want the charged 
Higgs contribution to reduce Rf M by more than la, we require M H ± > 60 (220) GeV for 
| A w | = 50 (70). 

Since Rl xp is only +la away from Rb M , it is not necessary to keep the narrow window 
of Ma and M h o if we allow 2a below the experimental data. In this case, M h o and Ma can 
be widened to much larger masses, and so the po constraint on the ceiling of the charged 
Higgs mass will also be relaxed. However, Mn± cannot be too small otherwise Rb will be 
decreased to an unacceptable value. 
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Summarizing this section the constraints by B° — B° mixing, p , and Rb give the following 
preferred scenario: 

1. Ma — M h o = 80 - 120 GeV; 

2. |A M | ~ 50; 

3. |A tt | <0.3; 

4. 80 GeV < M H ± < 200 GeV. 



VI. Conclusions 

We have demonstrated in model III of the general two-Higgs-doublet model the charged- 
Higgs-fermion couplings can be complex, even in the simplified case of no tree-level FCNC 
couplings. The phase angle in the complex charged-Higgs-fermion coupling determines the 
interference between the standard model amplitude and the charged-Higgs amplitude in the 
process of b — > sj. We found that for | A t i | ~ 1 there is a large range of the phase angle 
(9 50° — 90° and 270° — 310°) such that the rate of b — > sj is within the experimental 
value for all range of M#± . In other words, the previous tight constraints on M#± from the 
CLEO b — ► 57 rate is relaxed, depending on this phase angle. In addition, we also examined 
the effect of this phase angle on the neutron electric dipole moment and discussed other 
experimental constraints on model III. The necessary constraints are already listed at the 
end of the last section. Here we offer the following comments: 

1. The phase angle induces a CP- violating chromoelectric dipole moment of the 6-quark, 
which leads to a substantial enhancement in neutron electric dipole moment. The 
experimental upper limit on neutron electric dipole moment thus places a upper bound 
on the couplings: | A« A^l sin 6* < 0.8 for M H ± « 100 GeV. This bound has large 
uncertainties due to the hadronic matrix element of the neutron and the factorization 
scale. 

2. The phase angle will also cause other CP- violating effects in other processes, e.g., the 



decay rate difference between b — > S'j and b — > S7 [0], and in lepton asymmetries of 



b — ► si + t . These processes will soon be measured at the future B factories. 

3. Other experimental measurements, like F° — F° mixing, p , and R^, constrain only the 
magnitude of the couplings and the Higgs-boson masses but not the phase angle. 

4. The B° — B° mixing measurement can only constrain the charged-Higgs mass and 
\X U \ loosely because the mixing parameter x d depends on \V t d\, which is not yet well 
measured. Other uncertainties come from the hadronic factors: Bg, and t)b- 
Actually, the mixing parameter Xd is often used to determine \V t d\- 

5. As we have mentioned, if Rf M gets closer to the SM value the constraint on the neutral 
Higgs-boson masses: Ma and M h o = 80 — 120 GeV will go away completely. On the 
other hand, the charged-Higgs boson mass is still required to be larger than about 60 
GeV (for jA^I = 50) in order not to decrease Rb significantly. 
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Figure 1: Contour plot of the branching ratio b — > 57 versus and the phase of Xtt^bb for 
various values of |A#A&&| = 3, 1, 0.5. The shaded areas are excluded by the NEDM constraint 
\d n \ < 10" 25 e-cm. 
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Figure 2: Contour plot of the B° — B° mixing parameter x<i in the plane of \X a \ and the 
charged Higgs boson mass. The experimental value is = 0.734 ± 0.035. 
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M h °=M A (GeV) 



+ 1 o ~ R e *P =0.21656 




M H + (GeV) 

Figure 3: The Rb = T(Z — > b^/T^d due to the presence of (a) the neutral Higgs bosons, 
Ma — Myfi and (b) the charged Higgs boson. 
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